SHEP 93/94-15 



Bottom- Tau Yukawa Unification in the 
Next-to-Minimal Supersymmetric Standard Model. 

B. C. Allanach and S. F. King 

Physics Department, University of Southampton 
Southampton, S09 5NH, U.K. 



Abstract 



We discuss the unification of the bottom quark and tau lepton Yukawa couplings 
within the framework of the next-to-minimal supersymmetric standard model. We 
compare the allowed regions of the m 4 -tan (3 plane to those in the minimal supersym- 
metric standard model, and find that over much of the parameter space the deviation 
between the predictions of two models is small, and nearly always much less than the 
effect of current theoretical and experimental uncertainties in the bottom quark mass 
and the strong coupling constant. However over some regions of parameter space 
top-bottom Yukawa unification cannot be achieved. We also discuss the scaling of 
the light fermion masses and mixing angles, and show that to within current uncer- 
tainties the results of recent texture analyses performed for the minimal model also 
apply to the next-to-minimal model. 



It was realized some time ago that the simplest grand unified theories (GUTs) 
based on SU(5) predict the Yukawa couplings of the bottom quark and the tau lepton 
to be equal at the GUT scale 

h{M GUT ) = X T (M GUT ) (1) 

where M GUT ~ 10 16 GeV. Q Assuming the effective low energy theory below Mq UT 
to be that of the minimal supersymmetric standard model (MSSM) the boundary 
condition in Eq.p] leads to a physical bottom to tau mass ratio mb/m T in good agree- 
ment with experiment [@]. Spurred on by recent LEP data which is consistent with 
coupling constant unification, the relation in Eq.|I]has recently been the subject of in- 
tense scrutiny using increasingly sophisticated levels of approximation ||, £§]. These 
analyses showed that for given values of the strong coupling constant a 3 (M z ) and 
nib, there are only two allowed regions of tan/? for each choice of top quark mass m t : 
a small tan/3 branch and a large tan/3 branch. [] However the results are strongly 
dependent on a.3 (Mz) and nib, as well as GUT scale threshold effects. 

The recent investigations of the relation in Eq.|I| have focused on the MSSM. In 
this paper we shall instead assume that the effective low energy theory below the 
GUT scale is the next-to-minimal supersymmetric standard model (NMSSM) rather 
than the MSSM. The NMSSM || involves a single gauge singlet iV and has the [i 
parameter set to zero, with its effect being replaced by terms in the superpotential 
like XNH1H2 and kN 3 , where A and k are dimensionless couplings particular to the 
NMSSM. The scalar component of N is assumed to develop a weak scale VEV, x, 
whose value will not directly enter our calculations, since N does not directly couple 
to quarks and leptons. The NMSSM is an alternative to the MSSM which is equally 
consistent with coupling constant unification, and is at least as well motivated as 
the minimal model. Since the renormalisation group (RG) equations of the heavy 
fermions involve the new couplings A and k which may be quite large, there is no 
reason to expect the allowed regions of the m t — tan {3 plane, consistent with Eq.|T], 
to bear any resemblance to those in the MSSM. However, as it turns out the allowed 
regions in the NMSSM are quite similar to those in the MSSM although for large 
values of k(Msusy) and X(Msusy) the large tan/3 branch cannot be achieved. 

1 The relation in Eq.[l] is not exclusive to minimal SU(5), but also applies to other simple GUTs 
such as 50(10) and E 6 , providing the Higgs doublets are embedded in the smallest representations, 
as in the minimal SU (5) model. 

2 Recall that the MSSM is a two Higgs doublet model where H\ is the Higgs doublet which 
gives mass to down-type quarks and charged leptons, while Hi gives mass to up-type quarks. The 
superpotential contains a term [iHiHi, and the ratio of the two Higgs vacuum expectation values 
(VEVs) is tan/3 = v 2 /v x . 
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Our procedure closely follows that of ref.H. We preface our discussion with a 
brief summary of our approach and approximations. For a grid of tan/? and m< 
values we determine the Yukawa couplings of the heavy fermions at the scale Msusy- 
The three gauge couplings g>j (i = 1, 2, 3) at Msusy are determined by running them 
up from their measured values at Mz- Having chosen values of A and k at Msusy 
we run all the couplings up to Mqut = 10 16 GeV using the SUSY RG equations and 
check if Eq.|l] is satisfied to an accuracy of 0.1%, for all the values of tan/3 and m t 
in the grid. Since our goal is to compare the results of the NMSSM to those of the 
MSSM at the same level of approximation, it is sufficient to work to one loop order 
in the RG equations between Msusy and Mqut- Similarly, we take Msusy = m t for 
convenience, and ignore all low energy threshold effects i.e. assume all SUSY partners 
are degenerate with the top quark. However we shall consider the important effects 
of GUT scale threshold effects by considering 

Xb(M GUT ) = 0.9X T {M GUT ) (2) 

rather than Eq.|l]. In some models such as flipped SU(5) Eq.|2| may apply at tree- 
level. 

The various couplings at Msusy — m t were determined as follows. The run- 
ning masses of the fermions rrif(fi), where f — b,r and fi is the scale, were de- 
termined by running them up from their mass shell values mj(mj) with effective 
3 loop QCD (g) 1 loop QED 0, p. This enables the Yukawa couplings to be deter- 
mined at m t by: 

x / \ V2m t (m t ) ( . 

vsmp 

h(m t ) = — (4) 

r/bV cos p 

\ ( \ v / 2m T (m T ) 

At (mt) = —. (5) 

7] T V COS p 



where 



m f (m t ) 



and the VEV v — y v\ + v£ — 246 GeV. Clearly the Yukawa couplings at m t depend 
on both a^(Mz) and m^rrib). In the NMSSM the additional parameters X(m t ) and 
k(mt) are unconstrained, and may be regarded as additional free input parameters in 
our analysis. Finally the gauge couplings at Msusy = TOt were determined from some 
input values at M z , by using the standard model RG equations (including 5 quark 
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flavours and no scalar fields). The input values were taken to be a^MzY 1 = 58.89, 
a 2 (M z )~ 1 = 29.75, with a 3 (M z ) = 0.10 - 0.12. Note that whereas the m t referred to 
here is always the running one, it can be related as in B to the physical mass by 

4 
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(7) 



Given the dimensionless couplings at Msusy = mt, they are then run up to 
M GUT = 10 w GeV using the following SUSY RG equations relevant for the NMSSM, 
which we obtained in a straightforward way from ref. ||. Including the full Yukawa 
matrices we find 
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where U, D and E are the up, down and charged lepton Yukawa matrices respectively 
and the GUT normalisation convention of g\ has been used. Dropping small Yukawa 
couplings Eq.|] reduces to 
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(9) 



Our results are displayed in Fig.l as contours in the tan/3 — m t plane consistent 
with Eq.0. We take a 3 (M z ) = 0.11, m b = A.25GeV and the NMSSM parameters 
X(m t ) and k(m t ) as indicated. The MSSM contour is shown for comparison and 
is indistinguishable from the NMSSM contour with A (mt) = 0.1 and k(m t ) = 0.5. 
In fact our plot for the MSSM based on 1-loop RG equations is very similar to 
the 2- loop result in ref.H]. The deviation of the NMSSM contours from the MSSM 
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contour depends most sensitively on X(m t ) rather than k(mt). Two of the contours 
are shortened due to either A or k blowing up at the GUT scale. For A (m t ) = 0.5, 
k(m t ) = 0.5, no points in the m t — tan/3 plane are consistent with Eq|l| Yukawa 
unification, while for A (m t ) = 0.1, k (m t ) = 0.1 — 0.5 the contours are virtually 
indistinguishable from the MSSM contour. In general we find that for any of the 
current experimental limits on 0:3 and m&, the maximum value of X(m t ) or k(m t ) is 
~ 0.7 for a perturbative solution to Eq.[l|. 

In Figs. 2 and 3 we examine the effect of the experimental uncertainties in m& and 
a^[Mz) as well as the theoretical uncertainties parameterised by Eq.[|. In Fig. 2, for 
a rather low value of a^iMz) = 0.10, the region between the two solid lines respects 
Eq.[l| in the NMSSM with A(m t ) = 0.5, k = 0.1, and m b = 4.1 - AAGeV, while the 
dashed line satisfies Eq.Q for rrif, = AAGeV. The corresponding dashed line which 
satisfies Eq.^| for rrib = 4.1GeV occurs for m t < lOOGeV and so is to the left of 
the figure. Clearly for this value of strong coupling, the error in the bottom quark 
mass turns the sharp contours in Fig.l into thick bands, leading to no prediction of 
tan/3 for m t > 150GeV. In Fig. 3, for a higher value of a^iMz) = 0.12, the same 
uncertainties have a much less severe effect than in Fig.2, and larger values of the 
top quark mass are permitted. However the moral of Figs. 2 and 3 is clear: the sharp 
contours of Fig.l must not be taken too seriously given the present experimental and 
theoretical uncertainties. 

So far we have discussed the heavy third family fermion masses only. Let us now 
extend our discussion of the NMSSM to include the light fermion masses and mixing 
angles. It is well known in the MSSM that to one loop order in the RG equations, 
the running of the physically relevant Yukawa eigenvalues and mixing angles can be 
expressed in simple terms as shown below, 
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with identical scaling behaviour to V& of V u b, V ts , V tc i, where 

JlnM.su SY V 47r / 

To a consistent level of approximation V us , V u d, V cs , V c d, V t b, A u /A c , \d/K and A e /A M 
are RG invariant. The CP violating quantity J scales as V^ b . The Eqs. [10|, |TT1 
also apply to the NMSSM since the extra A and k paramters cancel out of the RG 
equations in a similar way to the gauge contributions as can easily be seen from Eq.^. 
The only difference to these physically relevant quantities is therefore contained in 
I T , I b and I t - 

In Figs. 4-6 we plot the values of I T , lb and I t as evaluated along the contours in 
Figs. 1-3 which satisfy Eq.|l|. Fig.4 illustrates the difference in the quantities between 
the MSSM and the NMSSM for the contour in Fig.l corresponding to X(m t ) = 0.5, 
k(mt) = 0.1. The NMSSM results are the upper lines of each pair, and it is clear 
that the deviation between the two models is small. In Fig. 5 we plot the three 
integrals for the NMSSM along the two solid contours of Fig. 2, corresponding to a 
low value of a^(Mz) = 0.10 and a range of nib = 4.1 — AAGeV . The experimental 
uncertainty in nib shown in Fig. 5 clearly swamps the theoretical difference between 
the values of and I t for the NMSSM and the MSSM shown in Fig.4, for this choice 
of as(Mz)- In Fig. 6 we plot the three integrals along the two solid contours of 
Fig. 3, corresponding to a higher value of a^iMz) = 0.12 and the same range of nib- 
Again the large deviations in lb and I t due to the error in nib swamp the theoretical 
differences between the two models. However, the values of I T in Figs. 5, 6 are quite 
robust, and the theoretical deviation in I T shown in Fig.4 is significant, and may play 
an important role in distinguishing between the two models. 

We emphasise that the results of the /j integrals shown in Figs. 4-6 play a key role 
in determining the entire fermion mass spectrum via the scaling relations shown in 
Eq .[To|. The small deviation between the NMSSM and the MSSM results compared to 
the experimental uncertainties, means that the recent GUT scale texture analyses of 
the quark mass matrices which were performed for the MSSM are equally applicable to 
the NMSSM. For example, the recent Ramond, Roberts and Ross (RRR) JTOf texture 
analysis is also based upon Eq.|l| and assumes a Georgi-Jarlskog (GJ) fli~i"| , |i~2| ansatz 
for the charged lepton Yukawa matrices, although their results in the quark sector are 
insensitive to the lepton sector. It is clear that all the RRR results are immediately 
applicable to the NMSSM without further ado since the only difference between the 
two models enters through the scaling integrals I{ whose deviation we have shown to 
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be negligible compared to the experimental errors. 

Finally in Fig. 7 we examine the question of full top-bottom-tau Yukawa unifi- 
cation in the NMSSM. Third family Yukawa unification has recently been studied 
in some detail in the MSSM 0] and occurs theoretically in minimal SO (10) |l4j 
and SU(4) <8> SU(2) 2 JL5J models. Fig.7 shows the values of the Yukawa couplings 
evaluated at the GUT scale for values of tan/3 along two of the contours in Fig.l cor- 
responding to the NMSSM with X(m t ) = 0.5 and k(m t ) = 0.1,0.4. The tau Yukawa 
coupling is of course equal to the bottom Yukawa coupling at the GUT scale and 
so is not labelled explicitly in this figure. In Fig.l it was observed that the high k 
lines are shorter than the others and consequently the large tan /3 regions cannot be 
achieved. The reason is clear from Fig.7, which shows that k(Mcur) blows up for 
large tan/3. This is because, as can be seen from Eq.[|, k scales quickly and readily 
becomes nonperturbative for k(m t ) > 0.5. The stunted lines in Fig.7 mean that, for 
k(mt) = 0.4 and the other parameters as assumed in the figure, top-bottom Yukawa 
unification cannot be achieved. The longer lines in Fig.7 show that, for k(m t ) = 0.1 
and all the other parameters unchanged, larger values of tan (3 may be achieved and 
top-bottom unification, which occurs for tan (3 ~ 50 for this set of parameters, is once 
again possible. 

In conclusion, we have discussed the unification of the bottom quark and tau 
lepton Yukawa couplings within the framework of the NMSSM. By comparing the 
allowed regions of the rat-tan j3 plane to those in the MSSM we find that over much 
of the parameter space the deviation between the predictions of two models which 
is controlled by the parameter A is small, and always much less than the effect of 
current theoretical and experimental uncertainties in the bottom quark mass and the 
strong coupling constant. We have also discussed the scaling of the light fermion 
masses and mixing angles, and shown that to within current uncertainties, the results 
of recent quark texture analyses []T0j performed for the minimal model also apply 
to the next-to-minimal model. There are however two distinguishing features of the 
NMSSM. Firstly, the scaling of the charged lepton masses will be somewhat different, 
depending on A and k. Although this will not affect the quark texture analysis of 



RRR, it may affect the success of the GJ ansatz fil| , fll^l for example. Secondly, the 
larger tan (3 regions may not be accessible in the NMSSM for large values of A and k, 
so that full Yukawa unification may not be possible in this case. 
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Figure Captions 



Figure 1. Contours in the m t — tan (5 plane over which the bottom-tau unification 
condition Eq.|l] is satisfied in the MSSM and the NMSSM. Central values of a 3 (Mz) = 
0.11 and m b = 4.25 GeV are assumed, and NMSSM contours for various X(m t ) and 
k(m t ) values are shown. 

Figure 2. Contours in the m t — tan/3 plane over which Eq.|I| is satisfied in the 
NMSSM with \{m t ) = 0.5 and k{m t ) = 0.1, for a low value of a 3 (M z ) = 0.10. The 
region between the two solid lines is for exact bottom-tau unification, with mj = 
4.1 - AAGeV. The dashed line satisfies X b = 0.9A r as in EqfJ for m b = AAGeV. 

Figure 3. Contours in the m t — tan f3 plane over which Eq.[I] is satisfied in the 
NMSSM with X(m t ) = 0.5 and k(m t ) = 0.1, for a high value of a 3 (M z ) = 0.12. 
The region between the two solid lines is for exact bottom-tau unification, with m b = 
4.1 — AAGeV. The region between the two dashed lines satisfies X b = 0.9A T as in Eq.||] 
for m b = 4.1 - 4.4GW. 

Figure 4. The If integrals defined in the text as evaluated along two of the 
bottom-tau unification contours in Fig.l. As in Fig.l, a 3 (Mz) = 0.11 and m b = 4.25 
GeV. The pairs of lines shown in this figure correspond to the MSSM (lower lines) 
and the NMSSM with \{m t ) = 0.5, k{m t ) = 0.1 (upper lines). 

Figure 5. The If integrals defined in the text as evaluated along the two exact 
bottom-tau unification contours of Fig.2. As in Fig. 2, a 3 {Mz) = 0.10, X b (m t ) = 0.5, 
k(mt) = 0.1. The shorter lines in this figure correspond to m b = A.lGeV, the longer 
lines to m b = AAGeV. 

Figure 6. The // integrals defined in the text as evaluated along the two exact 
bottom-tau unification contours of Fig.3. As in Fig. 3, a 3 (Mz) = 0.12, X b (mt) = 0.5, 
k(m t ) = 0.1. The shorter lines in this figure correspond to m b = AAGeV, the longer 
lines to m b = AAGeV . 

Figure 7. The various NMSSM couplings evaluated at the GUT scale, corre- 
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sponding to the two bottom-tau unification contours in Fig.l with X(m t ) = 0.5. The 
longer lines in this figure are from the k{m t ) = 0.1 line, and exhibit top-bottom uni- 
fication for tan (3 « 50. The shorter ones are from the k(m t ) = 0.4 curve, and show 
that top-bottom unification cannot be achieved since k blows up before a sufficiently 
large tan (3 can be achieved. 
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